JJSA0AC  TECHNICAL  LIBRARY 


5  0712  01020552  3 


MISCELLANEOUS  PAPER  N*74-l 


DECOUPLING  OF  GROUND  SHOCK  FROM 
EXPLOSIONS  IN  ROCK  CAVITIES 

by 

J.  L  Drake 


I 

I 


I 


0  i 


1 


January  1974 

Sponsored  by  Office,  Chief  of  Engineers,  U.  S.  Army 

Conducted  by  U.  S.  Army  Engineer  Waterways  Experiment  Station 
Weapons  Effects  Laboratory 
Vicksburg,  Mississippi 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originotor. 


The  findings  in  this  report  ore  not  to  be  construed  os  on  officiol 
Deportment  of  the  Army  position  unless  so  designoted 
by  other  outhorized  documents. 


MISCELLANEOUS  PAPER  N-74-1 


DECOUPLING  OF  GROUND  SHOCK  FROM 
EXPLOSIONS  IN  ROCK  CAVITIES 

by 

J.  L  Drake 


January  1974 

Sponsored  by  Office,  Chief  of  Engineers,  U.  S.  Army 
Project  A880 


Conducted  by  U.  S.  Army  Engineer  Waterways  Experiment  Station 
Weapons  Effects  Laboratory 
Vicksburg,  Mississippi 

ARMY-MRC  VICKSBURG.  MISS 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Foreword 


This  report  is  part  of  a  paper  presented  at  the  Department  of 
Defense  Explosives  Safety  Board  Sponsored  15th  Annual  Explosive  Safety 
Seminar  held  in  San  Francisco,  Calif. ,  on  18-20  Sep  1973-  The  reported 
research  is  a  generalization  of  a  formulation  of  a  model  for  stress 
wave  propagation  in  jointed  rock  which  was  sponsored  by  the  Office, 

Chief  of  Engineers  (OCE),  Military  Construction  Directorate,  under  Nu¬ 
clear  Construction  and  Engineering  Project  A880,  work  unit,  "Vulnera¬ 
bility  of  Deep  Underground  Structures  in  Rock,"  which  was  monitored  by 
Mr.  D.  S.  Reynolds,  OCE. 

This  report  was  written  by  Mr.  J.  L.  Drake  under  the  general  super¬ 
vision  of  Mr.  L.  F.  Ingram,  Chief,  Phenomenology  and  Effects  Division 
and  Mr.  W.  J.  Flathau,  Chief,  Weapons  Effects  Laboratory,  U.  S.  Army  En¬ 
gineer  Waterways  Experiment  Station  (WES). 

COL  G.  H.  Hilt,  CE,  was  Director  of  WES  during  the  preparation  and 
publication  of  this  report.  Mr.  F.  R.  Brown  was  Technical  Director. 
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Initial  cavity  radius 
Compress ional  wave  speed 
Detonation  velocity 
Arbitrary  function 
A  constant 
Detonation  pressure 
Radial  coordinate 
Dimensionless  radial  coordinate 
Time 

Radial  particle  displacement 
Radial  particle  velocity 
Peak  particle  velocity 

Peak  particle  velocity  of  cavity  wall  for  tamped  explosion 
Particle  velocity  of  cavity  wall 
Dimensionless  particle  velocity 
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Time  constant 
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Summary 


Decoupling  explosive  energy  by  detonating  charges  in  cavities 
larger  than  the  charge  is  an  effective  method  for  reducing  ground  shock 
magnitudes  over  those  of  fully  contained  bursts.  The  blasting  industry 
deliberately  decouples  explosives  in  presplitting,  smooth  wall,  and 
cushion  blasting  operations.  This  method  has  been  explored  as  a  pos¬ 
sible  method  for  concealment  of  underground  nuclear  explosions.  Decou¬ 
pling  can  be  an  effective  means  of  significantly  reducing  the  energy 
coupled  into  the  earth  by  accidental  explosion  of  magazines. 

This  paper  presents  an  analysis  of  ground  motions  generated  by  de¬ 
coupled  explosions  in  rock  cavities.  A  simplified  elastic  solution  is 
used  to  calculate  particle  motion  magnitudes  and  time  histories  produced 
by  an  explosion  in  rock  as  a  function  of  the  initial  cavity  radius  and 
the  loading  density  of  the  explosive.  Loading  density  is  defined  as  the 
total  explosive  weight  divided  by  the  initial  cavity  volume.  Initial 
cavity  conditions  were  estimated  from  the  adiabatic  expansion  of  the  ex¬ 
plosive  products  as  an  ideal  gas  for  slightly  decoupled  conditions,  and 
empirical  data  were  used  to  extrapolate  initial  cavity  conditions  to  low 
loading  densities. 

Results  compare  favorably  with  data  obtained  from  decoupling  ex¬ 
periments  in  three  different  rock  types  where  sources  ranged  from  small 
high-explosive  charges  to  low-yield  nuclear  explosions.  Spatial  attenu¬ 
ation  of  peak  particle  velocity  or  strain  and  positive  duration  of  par¬ 
ticle  motion  are  dependent  upon  the  initial  cavity  size  and  magnitude  of 
motion  on  the  loading  density. 
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DECOUPLING  OF  GROUND  SHOCK  FROM  EXPLOSIONS 


IN  ROCK  CAVITIES 


Introduction 


1.  Decoupling  explosive  energy  by  detonating  charges  in  cavities 
larger  than  the  charge  is  an  effective  method  for  reducing  ground  shock 
magnitudes  over  those  of  fully  contained  bursts.  The  blasting  industry 
deliberately  decouples  explosions  iin  presplitting,  smooth  vail,  and 
cushion  blasting  operations.  This  method  has  been  explored  as  a  pos¬ 
sible  method  for  concealment  of  underground  nuclear  explosions.  Decou¬ 
pling  can  be  an  effective  means  of  significantly  reducing  the  energy 
coupled  into  the  earth  by  accidental  explosions  of  magazines. 

2.  This  report  presents  an  analysis  of  ground  motions  produced  by 
decoupled  explosions  in  rock  cavities.  A  simplified  elastic  solution  is 
used  to  calculate  particle  motion  magnitudes  and  time  histories  produced 
by  an  explosion  in  rock  as  a  function  of  the  initial  cavity  radius  and 
the  loading  density  of  the  explosive.  Initial  cavity  conditions  were 
estimated  from  the  adiabatic  expansion  of  the  explosive  products,  and 
empirical  data  from  decoupling  experiments  were  used  to  extrapolate 
initial  cavity  conditions  to  lov  loading  densities.  This  presentation 
is  not  intended  to  be  a  rigorous  solution  of  the  decoupling  problem  but 
rather  a  semi-empirical  method  which  is  relatively  simple  and  accurate 
for  estimating  ground  motions  from  explosions  in  cavities. 


Elastic  Solution 


3.  The  general  solution  to  the  elastic  wave  equation  in  spherical 


coordinates  is  (see  reference  1  for  example): 


(i) 
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where 


u  =  radial  particle  displacement 
r  =  radial  coordinate 
t  =  time 

a  =  initial  cavity  radius 
c  =  compressional  wave  speed 

F(t)  is  an  arbitrary  function  of  the  argument  t  =  t  -  [jr  -  a)/cj  , 
which  is  the  reduced  time.  Expressions  for  particle  strain,  velocity, 
stress,  etc.,  can  be  written  directly  from  equation  1. 

k .  Particle  velocity  is  used  in  this  study  because  it  is  the  most 
commonly  measured  and  best  understood  ground  motion  parameter.  It  can 
be  determined  by  differentiation  of  equation  1  as: 


v(r,t)  =  - 


(2) 


where  a  primed  notation  indicates  differentiation  with  respect  to  the 
argument.  This  expression  is  an  ordinary  differential  equation  which 
can  be  easily  integrated  to  determine  F(x)  if  the  particle  motion  is 
specified  at  some  point  in  the  medium.  Assuming  the  particle  velocity 
at  the  cavity  boundary  is  v(a,t)  =  v^(t)  ,  then  integration  of  equa¬ 
tion  2  gives: 


F'(t) 


(3) 


Because  of  the  sensitivity  of  the  decoupling  phenomenon  to  the  initial 
cavity  radius,  it  is  convenient  to  integrate  equation  3  by  parts  to  ob¬ 
tain  a  series  expansion  in  powers  of  the  small  term  (a/c).  This  pro¬ 
cedure  greatly  simplifies  the  resulting  expressions  for  estimating 
particle  velocity  at  large  distances.  Noting  that 
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assuming  the  medium  was  initially  quiet  at  t  <_  0  ,  a  formula  for  par- 
tide  velocity  can  easily  he  found  hy  substituting  the  series  in  equa¬ 
tion  h  into  equations  2  and  3  as: 


v(r,t)  =  (r)  [v0(t)  +  ~  V0(t)  +  (5) 

5.  Previous  investigations  (reference  2)  have  shown  that  an  input 
particle  velocity  function  of  the  form 


to(t)  =  voe 
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was  adequate  for  describing  particle  motions  from  contained  explosions. 
Here  it  is  assumed  that  the  peak  particle  velocity  at  the  cavity  bound¬ 
ary  Vq  and  the  time  constant  a  are  functions  of  the  loading  density 
of  the  explosion.  Substituting  equation  6  into  equation  5  gives: 
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Note  that  the  peak  particle  velocity  occurs  at  x  =  0  ,  which  is  the 
arrival  of  the  wave  at  radius  r  and  is  given  by: 
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Initial  Cavity  Conditions 


6.  In  general,  initial  cavity  conditions,  v^  and  a  ,  are 
functions  of  the  loading  density  and  detonation  characteristics  of  the 
explosive  as  well  as  the  elastic  constants  of  the  rock.  The  initial 
cavity  velocity  is  proportional  to  the  peak  gas  pressure  in  the  cavity. 
Assuming  that  the  adiabatic  expansion  of  an  ideal  gas  can  be  used  to 
estimate  the  pressure  in  the  decoupled  cavity,  the  peak  velocity  at  the 
cavity  boundary  is  given  by: 
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where 

v  ,v  (p  )  =  particle  velocity  parameters  for  the  fully  coupled  and 
decoupled  charges,  respectively 

Pq,p^  =  explosive  density  and  loading  density,  respectively 
Y  =  ratio  of  specific  heats  of  the  explosion  gases 

Loading  density  is  defined  as  the  weight  of  the  explosive  divided  by 
the  cavity  volume. 

7.  The  particle  velocity  parameter  for  the  fully  coupled  condi¬ 
tion  can  be  estimated  from  the  detonation  pressure  and  the  characteris¬ 
tic  shock  wave  impedances  of  the  explosive  and  rock.  The  detonation 
wave  in  the  explosive  is  partially  reflected  and  partially  transmitted 
by  the  rock  medium.  Considering  this  reflection  process  to  be  acoustic, 
the  peak  pressure  a  in  the  rock  is: 
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detonation  pressure 

characteristic  impedances  of  the  rock  and  explosive, 
respectively 


The  impedance  of  the  explosive  is  taken  as  the  product  of  the  density 

p.  and  the  detonation  velocity  c  .  For  most  explosives,  the  detona- 
U  e 

tion  pressure  can  be  expressed  as: 


P 

e 


,  2 
kp0Ce 


(11) 


where  k  is  a  constant  less  than  unity.  Combining  equation  11  with 
equation  10,  v^  is  given  as: 
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(1 2) 


The  constant  k  was  determined  as 


k  =  |  (13) 

by  comparing  equation  8  with  experimental  peak  particle  velocity  data 
from  fully  coupled  experiments  in  granite  (reference  3),  salt  (refer¬ 
ences  h  and  5),  tuff  (reference  6),  and  alluvium  (reference  7)* 


Comparison  with  Decoupling  Experiments 


8.  It  is  convenient  to  introduce  the  dimensionless  scaled 
variables 


(14) 


and  rewrite  the  equation  for  peak  particle  velocity  (equation  8)  as: 
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This  expression  was  compared  with  data  from  decoupling  experiments  to 
establish  a  functional  form  for  a(p^)  and  to  verify  the  dependency  of 
the  particle  velocity  parameter  on  loading  density  as  given  by  equa¬ 
tion  9- 

9.  Relative  peak  strain  and  particle  velocity  data  from  experi¬ 
ments  in  granite  (reference  8),  limestone  (reference  8),  and  both  nu¬ 
clear  (reference  9)  and  conventional  (reference  10)  explosions  in  salt 
are  shown  in  fig.  1  as  a  function  of  the  relative  loading  density  of 
the  explosive.  These  data  follow  equation  9  quite  nicely  for  an  overall 
value  of  the  ratio  of  specific  heats  of  the  explosion  gases  of  1.2  over 
a  large  range  of  relative  density.  Considering  that  several  different 
explosive  types  (both  nuclear  and  chemical)  were  used,  the  scatter  in 
the  data  is  quite  acceptable. 

10.  Values  of  peak  particle  velocity  calculated  using  equation  15 
are  compared  with  experimental  data  from  the  Cowboy  (reference  10) 
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Fig.  1.  Relative  peak  particle  velocity  versus 
relative  loading  density 
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high-explosive  experiments  in  rock  salt  (see  fig.  2).  The  calculated 
curves  are  in  good  agreement  with  data  over  a  range  of  relative  loading 
densities  from  0.1  percent  to  fully  coupled.  Note  that  the  (a/r)  term 
becomes  more  and  more  dominant  as  the  relative  loading  density  decreases, 
indicating  the  increasing  influence  of  a  . 

11.  The  time  constant  parameter  a  is  difficult  to  determine 
with  any  degree  of  confidence.  Data  as  reported  by  the  Bureau  of  Mines 
(reference  8)  were  taken  from  long  cylindrical  explosives,  which  adds 
uncertainty  in  the  reported  periods.  Values  from  the  Cowboy  experiments 
(reference  U)  suggest  that 


a  = 


(16) 


which  is  verified  by  data  from  coupled  explosions  in  tuff  (reference  6). 


Conclusions 


12.  A  simplified  elastic  theory  is  presented  for  calculating  par¬ 
ticle  motion  magnitudes  and  time  histories  produced  by  an  explosion  in 
rock  as  a  function  of  the  initial  cavity  radius  and  loading  density  of 
the  explosive.  Initial  cavity  conditions  estimated  from  the  adiabatic 
expansion  of  the  explosive  products  compared  favorably  with  data  ob¬ 
tained  from  decoupling  experiments  in  three  rock  types  where  sources 
ranged  from  small  high-explosive  charges  to  low-yield  nuclear  explosions. 
An  average  ratio  of  specific  heats  of  the  explosion  gases  of  y  =  1.2 
was  found  to  fit  experimental  data  over  a  range  of  loading  densities 
from  0.05  to  100  percent  of  the  condensed  explosive.  The  scaled  time 
constant  was  found  to  vary  as  the  inverse  cube  root  of  the  relative 
loading  density. 

13.  Combining  the  results  of  the  elastic  calculation  for  peak 
particle  velocity  and  estimates  of  coupling  effects  on  initial  cavity 
parameters,  a  formula  for  predicting  particle  velocity  magnitudes  is 
found  as : 
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Fig.  2.  Peak  particle  velocity  for  Project  Cowboy 
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where  v^  can  be  calculated  from  the  detonation  velocity  and  initial 
density  of  the  explosive  by  the  expression  given  in  equation  12,  It  can 
be  seen  that  for  y  =  1.2  the  peak  particle  velocity  decreases  in¬ 
versely  as  the  cavity  radius  to  the  1.6  power  at  close  ranges  and  de¬ 
creases  inversely  as  the  cavity  radius  to  the  0.6  power  at  large  dis- 
tances  from  the  detonation.  The  l/r  term  becomes  more  dominant  for 
low  relative  values  of  the  loading  density. 

ik .  The  predicted  particle  velocity  time  history  has  an  abrupt 
rise  followed  by  an  exponentially  decaying  pulse  which  has  a  positive 
duration  of: 
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(18) 


which  has  a  component  due  to  the  cavity  size  and  a  second  term  which  de¬ 
creases  with  range  and  is  independent  of  the  initial  cavity  radius. 

This  result  is  verified  by  the  Cowboy  (reference  H)  and  Sterling  (refer¬ 
ence  9)  decoupling  experiments  in  salt  where  decoupled  periods  decreased 
to  the  value  of  a/c  .  A  rise  time  can  be  added  to  the  predicted  pulse 
without  difficulty.  It  is  showq,  in  reference  1  that  the  rise  time  pa¬ 
rameter  for  fully  coupled  charges  is  dependent  upon  the  rock  jointing 
patterns  and  is  independent  of  the  explosive  size. 
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tained  bursts.  The  blasting  industry  deliberately  decouples  explosives  in  presplit¬ 
ting,  smooth  wall,  and  cushion  blasting  operations.  This  method  has  been  explored  as 
a  possible  method  for  concealment  of  underground  nuclear  explosions.  Decoupling  can 
be  an  effective  means  of  significantly  reducing  the  energy  coupled  into  the  earth  by 
accidental  explosion  of  magazines.  This  paper  presents  an  analysis  of  ground  motions 
generated  by  decoupled  explosions  in  rock  cavities.  A  simplified  elastic  solution  is 
used  to  calculate  particle  motion  magnitudes  and  time  histories  produced  by  an  explo¬ 
sion  in  rock  as  a  function  of  the  initial  cavity  radius  and  the  loading  density  of 
the  explosive.  Loading  density  is  defined  as  the  total  explosive  weight  divided  by 
the  initial  cavity  volume.  Initial  cavity  conditions  were  estimated  from  the  adia¬ 
batic  expansion  of  the  explosive  products  as  an  ideal  gas  for  slightly  decoupled  con¬ 
ditions,  and  empirical  data  were  used  to  extrapolate  initial  cavity  conditions  to  low 
loading  densities.  Results  compare  favorably  with  data  obtained  from  decoupling  ex¬ 
periments  in  three  different  rock  types  where  sources  ranged  from  small  high- 
explosive  charges  to  low-yield  nuclear  explosions.  Spatial  attenuation  of  peak  par¬ 
ticle  velocity  or  strain  and  positive  duration  of  particle  motion  are  dependent  upon 
the  initial  cavity  size  and  magnitude  of  motion  on  the  loading  density. 
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